Spin stiffness calculation in anisotropic XY model with Ring exchange interaction. 
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We present the spin wave theory of XY model with anisotropic nearest neighbour (NN) interac- 
tions J(J') along the x(y) directions, next nearest (NNN) neighbour interaction Jd and the ring 
exchange interaction K on the square lattice. We calculate the thermodynamic quantities: Zero 
temperature spin stiffness, internal energy, specific heat and the magnetization. Using the diago- 
nalized Hamiltonian, we show that no soft modes develop when r\ + A > 0, where r\ — J' /J and 
A = K/J. We further show that anisotropy (rj — 2) decreases the spin stiffness by 5.7% of its 
isotropic (77 = 1) maximum value for some values of A and 8 — Jo/ 'J- A similar reduction shows up 
in the magnetization. The plot of the stiffness against r?(A) reaches a maximum at r\ = 0(A = 0) for 
specific values of S and decreases rapidly as it approaches 77 + A + 1 = 0. In general, the supersolid 
phase transition suggested earlier 10 will occur in the regime r\ + A + 1 < 0. 



I. INTRODUCTION 

The quantum XY model has been a subject of interest 
for many decades. It has be en intensively studied both 
numerically and theoreticall} |3 | 4 | 8 | 10 | 13 1M ^l. This model 
is of considerable importance for the understanding of 
quantum phase transitions and thermodynamic proper- 
ties at low temperatures. Most theoretical studies of this 
model is based on spin wave theory method. 

This method involves the mapping of the spin opera- 
tors to bosonic operators such that the spin commutation 
relation is satisfied. It provides a good approximation 
for obtaining the low-lying excitation spectrum of quan- 
tu m spin systems. Several versions of spin wave theory 
exit! 1 * 2 * 12 ! The most popular one is based on the Holstein- 
Primakoff representation^ which was first applied to the 
study of Heisenberg model by AndersorP and further ex- 
tended to second order by KubcP and OguchP^. 

Gomez-Santos and Joannopoulous^ showed how this 
method can be applied to the case of isotropic XY model 
by making a good choice of the quantization axis. Since 
then numerous applications of spin wave theory have 
been carried out on the isotropic XY model with differ- 
ent lattice configurations. Most of the results obtained so 
far are in a good agreem ent with quantum Monte Carlo 
simulations(QMC) ^ 13 ' 14 '. 

The isotropic XY can be supplemented with other 
non-trivial interactions such as anisotropic nearest neigh- 
bour (NN) interaction, next nearest neighbour (NNN) 
interaction or the ring exchange interaction. These in- 
teractions can lead to exotic quantum phases in the sys- 
tem such as superfluid* ^, valance-bond-solid (VBS) 25 
and gapless Mott insulator (GMI)P^ phases. They can 
as well destroy the phases depending the strength of the 
interaction. 

The ring exchange interaction was first introduced to 
study the magnetic properties of solid 3 He^. This ex- 
change interaction, alone or in competition with the pure 
XY model (NN exchange only) has been studied ex- 
tcnsivel y on th e square lattice using different numerical 
methodsPEsEU. A linear spin wave theory approxima- 



tion to this model has also been carried on the square 
latticeP. 

It is well known that the ring exchange interaction 
destroys the superfluid phase of isotropic XY modeP^. 
However, in this paper we shall consider the XY model 
with all the interactions mentioned above: Anisotropic 
NN, NNN and the ring exchange interactions. This in- 
troduces four (or three dimensionless ) parameters into 
the Hamiltonian in contrast to just one parameter. The 
format of the paper is as follows: In Sec. II, we present 
the model Hamiltonian and the classical ground state. 
In Sec. Ill, we apply spin wave theory by choosing our 
quantization axis along the x-direction. Diagonalize the 
Hamiltonian by means of Bogoluibov transformation and 
obtain the spectrum. In Sec. IV, we plot the spectrum 
and the magnetization as a function of the parameters 
involved. We show that the soft mode of the spectrum 
imposes a constraint on the parameters of the Hamilto- 
nian. In Sec.V, we calculate the zero temperature spin 
stiffness by applying a twist along the x-direction and 
plot it as a function of the parameters . Finally, in Sec. VI 
we make some concluding remarks. 



II. MODEL 

The model Hamiltonian we shall consider in the paper 
involves the four different interactions. The exchange in- 
teractions J(J') with summation over the NN along the 
x(y) directions, exchange interaction Jr> with the sum- 
mation over NNN along the diagonal and finally the ring 
exchange interaction K with summation over the square 
plaquettes. 

H = J E ( S t S F + M - J ' E ( S t S k + h.c) 

{ij) (J'fe) 

- Jn ]T (S+S- + h.c) (S+STS+S- + h.c) 

[ik] (ijki) 

(1) 
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The spin operators obey the usual commutation relation theory we have 



S?,S 
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± iff: 



iSa^jS^Sij, where a, /3, 7 = x, y, z and Sf = 

; . The model Hamiltonian has a sign problem 
in quantum Monte Carlo simulations in the regime K < 
0. However, this sign problem cannot be captured by a 
simple spin wave theory calculation. In the limit J = J' 
and = 0, the Hamiltonian reduces to the one studied 
in Ref. 10. When J, J' » Jd,K, this model undergoes 
a Kosterlitz-Thouless phase transition^ at Tkt ~ 0.69 
for 2D model and a superfluid phase for temperatures less 
thanT TK ESl. The regime Jjj — 0, K < has been studied 
in Ref. 23 using variational Monte Carlo (VMC) and 
density matrix renormalization group (DMRG) method. 
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FIG. 1: (Color online): Actions of the exchange interactions. 
The hoping strength J acts along the a;-axis, J' acts along 
the y-axis, Jd acts along the diagonal and K is the four spin 
cyclic ring exchange interaction. 



Ferromagnetic interaction corresponds to J,J',Jd,K > 
0. All calculations and plots will be done as a function 
of the dimensionless quantities: 77 = J' /J, 8 = Jd/J, 
and A = K/J. We will set J = 1/2 anywhere else which 
corresponds to the parameter value of pure XY model . 
The classical ground state is found by treating the spins 
as classical vectors i.e Sf = Se iS . For spin- 1/2 we have 



^ = -^(l + 77 + 2* + A/4). 



(2) 



where we have taken the coordinate number z — 2 for 
NN interactions and z = 4 for NNN interaction . 



Sf « X - ( a] - a, 



(3) 



The spatial Fourier transformation of the bosonic opera- 
tors is written aiP 

k k 

where the momentum k runs over the first Brillouin 
zone (BZ) of a square lattice with unit nearest-neighbour 
distance i.e n < k x < tt, tt < k y < ir. Inserting 
Sf = S$ ± iSy, ^ and Q into 0, the quadratic (non- 
interacting) part can be grouped in the following way 



Hn = H 



MF 



^ [-4k (a^ab + ai k a_ k 



-B k ( a k aL k + a k a_ k 



(5) 



This Hamiltonian is diagonalized by the Bogoluibov 
canonical transformation to quasiparticle boson opera- 
tors «k and ajp 112 * 

ak = IkUk - rn k al k , a k = £ k a k - m k a_ k , (6) 



with 



m k = 



A - £k 



(7) 



+ £k 

2e k ' '"" V 2e k 

Applying the above transformations, the diagonalized 
quadratic (non-interacting) part yields 

H = H M f + (£k - -4k) + £k(a k a k + aL k a_ k )- 

k k 

(8) 

The mean-field energy and the coefficients are given by 
-JN 



Hmf — 



(1 + rj + 26 + A/4) 



Ak = J [Qk + Ai? k ] , 
Bk = J [5k + AT k ] , 



£k 



(9) 
(10) 
(11) 

(12) 



III. SPIN WAVE THEORY 

The basic assumption of spin wave theory lies on se- 
lecting a classical ground state and determining the fluc- 
tuation around it. In other words, one considers quan- 
tum fluctuations very close to an ordered ground state 
configuration of the system under study. By choosing 
the quantization axis along the x-direction (instead of 
the z-direction), one can then write the spin operators in 
terms of the bosonic operators using the famous Holstein- 
Primakoff transformation^^]. In the linear spin wave 



where 



Sk 
Tk 



(l-^)+<l-^)+^(2-7^J 

7k, , 7fc 



1. 1 1 

-7ViK +7fcj,) - -ilk x lk v , R k =- 



and the lattice structure constants are 

j kx = cosfcz, 7fc H = cosfcy. 



(13) 
(14) 

Tk, (15) 
(16) 
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IV. THERMODYNAMIC PARAMETERS 

The diagonalized quadratic Hamiltonian ^ gives the 
spin wave ground state and the excitation spectrum. 
The energy spectrum is given by £k = 2et. Figure ^ 
shows the plot of the spin wave energy long the direc- 
tion k x = k y . In the top figure, the spectrum shows 
two zero modes at k = (0,0) and k = (w, n) for 77 = 1, 
A = — 2 and several values of 8. This result agrees with 
the result found previously, that no soft modes develop 
in the regime A > (for the isotropic case rj = 1 and no 
NNN interaction Thus, the NNN interaction does not 
change the soft modes of the energy spectrum but only 
increases its peak. This is because the ring exchange 
term has already incorporated NNN sites. The bottom 
figure is interesting, the zero (soft) mode of the energy 
at k = (tt, tt) is gapped as one moves away from rj = 1. 





FIG. 2: (Color online): The energy spectrum along k x = k y 
for rj = 1, A = —2 and several values of <5 (Top) and for several 
values of 77, A and S (Bottom). The destroyed soft mode at 
k = (tt, tt) for 77 > 1 is restored whenever 77 + A + 1 = 0. 



However, we found that this soft mode is restored 
whenever 77 + A = — 1 . To show how this condition came 
about, first consider the expansion of the energy spec- 
trum near the zero mode k = (0, 0): 



e(k -> 0) = [v 2 x kl 
The sound speeds v x and v y are expressed as 

v x = JV(2 + 45 + 277 + \)(28 + 1) 
v y = JV (2 + iS + 277 + A)(2£ + 77) 



(17) 

(18) 
(19) 



Thus, the spectrum is gapless at k = independent of the 
values of 77, A, and 8 as shown in Fig.pl). For a specific 
case 77 = 1, 8 = 0, and A = —2, we have 



e(k^ 0) = c|k| 



(20) 



where c = J\[2, Now consider the Taylor expansion of 
the energy spectrum near the zero mode k = (7r,7r): 



e(q -> 0) = [v 2 x q 2 x + v 2 y q 2 y + As] 



(21) 



where q x — k x — tt, etc. and the sound speeds v x and v y 
in this case are 



J 



V(2 + 45 + 2r7 + A)(2<5- A + 1) 



«„ = - V( 2 + 45 + 277 + A)(2<5 - 77 - A) 
The energy gap Ae is of the form 

Ae = J 2 (l + 77 + A)(2 + 48 + 2r] + A) 



(22) 
(23) 

(24) 



It's clear that at q = 0, no soft modes develop when 
77 + A > 0. This is illustrated in Fig. ^ (bottom). Notice 
that the spectrum is imaginary in the regime 77+A+l < 0, 
indicating a change in the ground state configuration. 

We will now consider the thermodynamic quantities: 
Total internal energy and specific heat capacity of the 
spin waves. At low temperature the spin wave are in- 
dependent bosons and only low energy spin waves are 
excited. We can then use Eqs.([T7l) and (f2ll) to determine 



the power of the temperature dependence on the specific 
heat C s . In the regime where Ae = 0, i.e 1 + 77 + A = 0. 
Eq. (21 ) becomes 



e(q) = [e 2 x + e 2 y ] 



2] 1/2 

y\ 



(25) 



where e x = v x q x etc. The total energy of thermally ex- 
cited spin waves is 



U = 



g(q)^q 

e e(q)/T _ 

1 

(2ir) 2 v x Vy 



1 

OO 



(26) 



de-r 



1 



The above integrals can be converted into polar coordi- 
nate which gives 



U 



1 



de- 



2ttv x v v J e £ / T - 1 
The spin wave specific heat yields 



C(3) 



T 3 

TTV x V y 



C, = 



dU 
df 



3C(3) T 2_ 



(27) 



Another thermodynamic quantity of interest is the mag- 
netization. Physical observables are calculated by intro- 
ducing the Green's function for the system. This is useful 
when one considers higher order expansions in 1 / S which 
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involve bosonic interactions. At zero temperature, the 
Green's function of the spin wave operator in this case is 



G U (M) = -t<Ta k (t)4(0)>, 



(28) 



where the subscript on the right hand side refers to the 
first diagonal term in the total Green's function 2 x2 
matrbP-H Using ([6]), the Fourier transform of (28) gives 

gum= , -V-Iri" v (29) 

The average magnetization is then expressed as 
1 / t \ 
1 1 ^ f°° rfw 



M 



(30) 



N ^ 



2iri 



Gn(k,w)e 



-iojQ 



Evaluating the second term by contour integration, we 
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FIG. 3: (Color online): Magnetization vs A for the isotropic 
case rj = 1 with S — O(black), S — 0.029(green) and the 
anisotropic case rj = 2 with S — O(blue), S — 0.029(red). The 
value of the magnetization for rj = 1, S = and A = is 
M = 0.439 and M = 0.434 for rj = 2, S and A remains the 
same. The decrease is 1.1% of M — 0.439. The magnetization 
drops rapidly when rj + A + 1 = 0. 



obtain the deviation 



AM 



L E 



2N 



Ah 



Bi 



(31) 



where and B^ are given by Eq.(10) and Eq.(lll re- 
spectively. In Fig. p]), the magnetization is plotted 
against A for the isotropic r\ = 1 and the anisotropic 
cases -q — 2 with two different values of 6. The mag- 
netization decreases rapidly when rj + A = —1 reflecting 
the fact that no soft modes is developed in the regime 
rj + A > 0. For <5 = and A = 0, anisotropy decr eases 
the value of the isotropic magnetization M = 0.43£P^MI 
by 1.1%. 



V. SPIN STIFFNESS 

The spin stiffness, or helicity modulus is the change in 
the ground state of a spin system resulting from apply- 
ing a twist <fi between every pair of neighbouring lattice 



bonds. It is given by the second derivative of the twisted 
ground state energy with respect to the twist 



Ps = 



1 d 2 E (fa 



N dfa 



(32) 



'»=() 



where E (fa is the ground state energy of the twisted 
Hamiltonian and N is the number of lattice sites of the 
system. In the thermodynamic limit, the sign of the spin 
stiffness determines if the system has a long-range mag- 
netic order(LRMO), p > means the existence of LRMO 
in the system whereas p s = means no LRMCB3. The 
above mathematical expression for p s is equivalent to the 
difference in the ground state energies 



E {fa-E 
N 



1 



(33) 



where E (fa = (H(fa) and E = (Ho). The twist depen- 
dent Hamiltonian is found by rotating the spin at site i 
by fa around the z-axis i.e Sf -> S+e 1 ^, Sr -> S^e'^, 
so H(fa) is given by 



flu>\ = - •/ > ] (S+Sre***-^ + /i.c) 
(S/Sfe e*to - **> +h.c\ 



(ij) 

r E 

(jk) 



Jd 



K 



[ik] 

E 

(ijki) 



(34) 



-<t>i+<i>k-<i>i) _|_ fa 



Now consider a uniform twist fa, along the positive x- 
axis nearest neighbour bonds. Using the labelling in 
fig.Q we have fa - fa = fa — fa - fa, 4>j - fa — 
= fa — fa, and fa — fa = <fi x . Thus, the twist de- 



pendence on the second and the fourth terms in (34) 



vanishes. MacLaurin expansion of (34) around <fi x = 
gives 



H(fa) = H 



1 



b x3 x 



?T S . 



(35) 



The spin current operator (paramagnetic term) and 
the spin kinetic energy operator (diamagnetic term) 
are defined mathematically as 



dH( 



x d(fi x 
Explicitly we have 



rpS 



d 2 H(fa) 



l+x 



(36) 



^l~^l+x+y 



h.c 



(37) 



i 



l+x—y 



h.c 



(38) 
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Using second order perturbation theory we obtairPSEI] 

(39) 



2 

N 



where |0) and \v) are the ground state and the excited 
states respectively. 



r 



at A = 0. This is obvious from the first term in (40 1. 
The anisotropy case rj = 2 is quite interesting. The spin 
stiffness decreases by 5.7% of its isotropic maximum value 
independent of the value of 6. The maximum values of 
p s still center at A = in this case. Notice that the 
spin stiffness shows a sharp decrease when rj + A + 1 = 
which indicates an ordered phase with (vr, tt) symmetry. 
The nature of this ordered phase for the case 77 = 1 and 
A = — 2 has been suggested to be a supersolid phased 
In general this phase transition will occur in the regime 
where 77 + A + 1 < 0. 

On the other hand, the plot of spin stiffness against 77 
shows a similar trend. It has a maximum at 77 = and 
decreases as 77 moves away from zero for some values of 
5. It also decreases so fast as it approaches the region 
where 77 + A + 1 = (see Fig.j5]|). 



V 

II 
II 



FIG. 4: (Color online): Spin stiffness vs A for isotopic case 
(top) n = l,6 = O(black), O.Ol(green), 0.015(blue), 0.02(red). 
Anisotropic case (bottom) 77 = 2, 8 the same. Anisotropy 
decreases the value of the spin stiffness. The stiffness shows a 
rapid decrease when the constraint 77 + A + l = 0is satisfied. 



At zero temperature the second term in (39) vanishes. 



The first term can be written in terms of the bosonic 
representations ^ and the Bogoluibov transformations 
(|6) which gives 



where the new coefficients arc 



(40) 



Ck = JQk(v = 0), D lc = JS K (r ! = Q). (41) 

Figure Q shows the plot of the spin stiffness against the 
parameter of the ring exchange A for two cases: Isotropic 
case (top) 77 = 1 and anisotropic case 77 = 2 (bottom) for 
different values of the NNN exchange 8. In the isotopic 
case we recover exactly the result obtained in Ref. 10 
for 77 = 1 and 5 = 0. As S increases above zero, the 
value of the stiffness increases with its maximum centers 



FIG. 5: (Color online): Spin stiffness vs 77 for A = 2, and 
S = 2(black), 2.01(green), 2.02 (blue), 2.03(red). The stiffness 
shows a rapid decrease when the constraint 77 + A + l = 0is 
satisfied. 



VI. CONCLUSION 

In this paper we have studied the hard-core boson 
( zero field XY model) using linear spin wave theory. 
In contrast to the previous study on this model, we in- 
cluded the effects of NN anisotropy, NNN and ring ex- 
change interactions which introduced three dimensionless 
parameter as opposed to just one. It was shown that 
no soft modes (Goldstone modes) develop in the regime 
77 + A > 0. The spin stiffness was obtained by apply- 
ing a twist along the x-axis nearest neighbour sites. We 
showed that anisotropy decreases the values of the stiff- 
ness by 5.7% of its isotropic maximum value for specific 
values of A and 5. Similar decrease was calculated for 
the magnetization. The stiffness shows a sharp decreases 
as it approaches 77 + A = — 1 which reflects the fact that 
no soft modes develops for 77 + A > 0. In general the 
supersolid phase with (77, tt) symmetry suggested in the 
previous work on this model will occur in the regime 
77 + A + l < 0. 
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